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PREFACE 


This report was prepared under study Contract 
NAS5-1 14-80, '‘Design Study of a Low Bit Rate nans - 
miner and Simulator, " by Scientific Management . 
Systems, Inc. The study was sponsored by the National 
Aeronautics and Space Administration under the direction 
of Richard L. Kutz. 

The purpose of this study was to establish, a oaseline 
communications channel model to stuay low bin rc.e 
telemetry transmissions from a Venus Planetary Explore! 
Spacecraft. The approach taken was to impiemeni 
existing models drawn from studies conducted by ihe 
Jet Propulsion Laboratories (JPL) and once these 
existing baseline simulation models were implemented 
and running on NASA/GSFC compuiers new and 
oriaina! concepts were to be aaded onto the baseline 
model as the Planetary Explorer Program requirements 
dictated. Subsequent to the completion of this study 
the Planetary Explorer Program was transferred >o 
Ames Research Center and the urgent requirements 
of this study have diminished. Consequently this 
report presents the results achieved prior to the 
termination of the Planetary Explorer Project at 
NASA/GSFC. 
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INTRODUCTION 


1.1 -GENERAL 

This report describes the results achieved to date in developing a simulation model of 
the Planetary Explorer Small Probe Communications Links. The results described are 
the first step in developing an overall model to be used as a test for evaluation of design 
alternatives for the PE communications subsystems. 


1.2 SIMULATION MODEL 

1.2.1 General 

The construction of appropriate mathematical models for physical sources and channels 
must be guided by the intended use of the finished product. It was decided to initially 
develop a model based upon low data rate signals which are recovered by spectral 
analysis using FFT techniques with signal acquisition and tracking being accomplished 
using a method implemented by JPL. Subsequent modifications and additions to the 
baseline simulation will be accomplished on an as required basis as dictated by required 
design tradeoffs. 

Figure 1 .1 presents a block diagram of a communications system slanted toward information 
theoretic considerations . ]/ These considerations when correlated with hardware consid- 
erations will provide a basis for detailed hardware design. The scope of this initial effort 
is limited to those blocks included in the dashed lines. Subsequent paragraphs describe 
the overall simulation concept with details for the discrete channel presented in section 2. 


1/ Information Theory & Reliable Communication , by R. G. Gal lager, 1968, 
John Wiley & Sons. 
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DISCRETE CHANNEL 
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Figure 1.1 BLOCK DIAGRAM OF COMMUNICATIONS SYSTEM 











1.2.2 Overall Communications System Channel Model 

The class of channel models which we are considering is the class where the set of inputs 
and outputs are each a set of time functions (waveforms) and for each input waveform the 
output waveform is a random process. For the initial application, the source. Figure 1.1 
is the small probe sensor complement and the source encoder is an A/D converter. When 
transmitting binary data over a channel in the above class it is convenient to separate the 
channel encoder and decoder each into two parts. The digital data coming out of the source 
encoder (A/D converter) enters the discrete channel encoder (for the small probe application 
a rate 5 convolutional encoder) producing a sequence of symbols from a fixed alphabet, 
say a ,...., a^at fixed time intervals, say one symbol each ^ seconds. The digital data 
(for our initial application an MFSK system) produces one of a, fixed set of 
waveforms, say S} (t), S»(t), . . . .S^(t), each of duration r^. Thus the entire waveform 
input to the channel has The form 


E 



(t-HT c ) 


0 ) 


where the sequence i , n = . . .- 1 , 0 , 1 , . . . is determined by the corresponding inputs to the ! 

digital data modulator. J 

Now this resultant waveform which is input to the channel is attenuated, modified in phase, j 

and Gaussian Noise is added as it passes through the channel. The digital data demodulator { 

takes the received waveform and converts it back to a sequence of symbols from a fixed alphabet. 

The discrete channel decodes and makes a decision as to what alphabet symbol has been trans- ■! 

• mitted in the corresponding time interval. In more sophisticated cases the output will also j 

contain information as to how rel iable the decision is, yielding an output alphabet of b^ ,...., b. . 

The source decoder and destination are functions of the data itself and will remain unspecified * 

for the purposes of this report. j 

j 

Since we have framed the overall models in an information theoretic sense we will be able | 

to make use of some standard measures of efficiency for these cases, i.e., channel capacity, j 

noisy channel coding theorem and desired bit rate. Subsequent paragraphs will describe j 

the modelling accomplished to date for the discrete channel portion of the overall model ! 

and the utility of these models. Future work will expand upon that discussed here with I 

emphasis upon the portions of Figure 2.2 not yet addressed and upon the output measures of I 

effectiveness. j 
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DISCRETE CHANNEL MODEL 


2.1 GENERAL 

The mathematical model for the discrete channel which has been developed is partitioned 
into two major groups, i.e., the input signal generation and the received signal acquisition. 
The approach taken has been that of establishing a baseline system utilizing the JPL design 
studies and the GSFC Phase A Explorer Summary Report. This baseline system is being 
used to Insure that the simulation concept is sufficiently flexible to permit detailed 
tradeoff analyses and that the concept incorporates provisions for growth, as required. 

2.2 SIMULATION LOGIC 

The gross logic of the simulation as presently envisioned is as follows. Steps a through d 
are contained in the input signal generator and e through g in the signal acquisition 
model. 

Signal Generation 

a. Establish transmitter parameters 

b. Compute coded MFSK data set 

c. Characterize signal at the receiver as a function of: 

(i) Path loss 

(ii) Random frequency instability 

(iii) Channel fading 

(iv) Channel phase variations 

d. Generate random variations on transmitted data set 

(Gaussian Noise) 
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Signal Acquisition 


e. Generate Fourier coefficients 

f. Establish frequency acquisition and synchronization 

g. Establish time acquisition and synchronization 

Thus, a baseline simulation which has been developed permits the generation of a 
simulated received signal which incorporates the effects of the transmitter parameters 
and channel parameters as identified in Table 3.1 . The effects of variations in these 
parameters are assessed by Monte Carlo methods. Increased sophistication in channel 
models can be accommodated without major modification to the baseline simulation 
which has been developed. 

The signal acquisition model is based upon the work previously accomplished at JPL. 

It permits the testing of such typical items as decision thresholds for the FFT outputs 
and acquisition times under various assumptions concerning the input signals and accuracy 
of the coded signal set. 

2.3 GROWTH ITEMS 

Growth items for the present baseline simulation include: 

a. Encoding and Decoding algorithms; 

b. Improved channel characterization; 

c. Alternate signal synchronization schemes; and 

d. Data window weighting of received signals. 
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INPUT SIGNAL GENERATION 


3.1 GENERAL 

The signal which is received at the digital data demodulator can be characterized as a 
function of the key transmitter parameters, the random disturbances in the channel and 
the effective noise temperature of the receiving system. The discussion which follows 
will incorporate into mathematical models the effects of the transmitter and channel 
key design parameters shown in Table 3.1 

3.2 INITIAL COMPUTATIONS 

The initial computations which will be made in the simulation are those of received 


( 2 ) 


(3) 

Pj, = Transmitter Power (watts) 

G ■ = Gain of the Transmitter Antenna (dimensionless) 


power and nominal signal to noise ratio as follows.!' 


P R= P , G , 


(4.8x1 0 11 )Lf 2 R 5 


S/N = PGD' 


(4.9)R T W 


where - Nominal Received Signal Power (watts) 


2/ Low Datg Rate Telemetry , Goldstein & Kendal, Proc. of Am. Astron. Soc. Symp. 
on Unmanned Exploration of the Solar System, February 1965. 
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TABLE 3.1 


KEY DESIGN PARAMETERS 



Number of Tones 
Tone Spectral Width 
(Osc. Thermal noise) 
Time Duration of Tones 
Coding Method 
Transmitted Power 
Tone Frequency Spacing 
Long Term Drift (Doppler 
residuals, long term 
changes, power supply 
variation) 

Carrier Frequency 
Bit Rate 
Antenna Gain 
Bits/Word 
Losses 


Planetary Motion (Venus) 
Wind Speed & Turbulence 
Doppler Shift 
Atmospheric Effects 
Gaussian Noise 
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Receiver 


Receiver Antenna Gain 
Receiver Noise Temperature 
Signal Detection Method 
Frequency Synchronization 
Method 

Time Synchronization Method 
Long term Drift Compensation 
Method 

Doppler Correction 
Oscillator Instability 
Decoding Method 









L = Losses in the Total System (dimensionless) 
f = Frequency of the Transmitter Carrier (hertz) 

R = Range from Transmitter to Receiver (Astronomical Units) 

S/N= Nominal Signal to Noise Power Ratio at the Receiver (dimensionless) 
D = Diameter of the Receiving Antenna (meters) 

T = System Effective Noise Temperature (degrees kelvin) 

W = System Bandwidth (hertz) 


3.3 INTRODUCTION OF RANDOM VARIATIONS 

Having established nominal values for transmitted power and signal/noise power ratio 
at the receiver, it is necessary to introduce the effects of random variations in the 
transmitter and channel parameters. After carrier demodulation the received signal 
will be assumed to be of the form: 


X(t) = A cos (2 7T ft-J$) +n(t) 


where: 


= Amplitude of the Received Signal 
= the signalling frequency plus frequency error 
= unknown phase term 
= random noise process 


In order to model the effects of the random variations we will use a Monte Carlo process 
in which probability distributions will be established for the random variables of interest. 
Random samples will then be drawn from th^se probability distributions for each of the 
random variables at times coinciding with the sampling rate of the system. These random 
samples will be added to the nominal values of the received signal in order to establish 
stochastic variations in the signal at each sampling interval. 


Mathematical !y, 

X(t.) = S(t.) + N(t.) 

S( t |> = j2(P R ^P ) cos 


N(t.) = 


2 < p r +ap r.> 

I 


1*2 JT f t.+ ^ 0 . + <& 9 . 1 
L K ' r i J 


cos ( A ) 
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where: 


S(t.) = 

N(t.) = 

AP R. = 

i 

f K = 
* 0 . = 
A9. = 

i 

A a = 


Received Signal at Time t. 

Noise Signal at Time t. 

Random Sample of Signal Amplitude Variation at Time t. 

Signalling Frequency 

Random Sample of Phase Variation due to Long Term Drift at Time t. 
Random Sample of Phase Variation Due to Short Term Drift at Time t. 
Random Sample of Phase Variation Due to Channel Noise at time t. 


Thus utilizing equations (2), (3), (5), (6) and (7) we have modelled the received signal in terms 
of the key design parameters for the transmitter and the channel shown in Table 3. 1 . Figure 3. 1 
presents a flow chart of the method by which the composite input signal for a specified time 
period will be computed in the software routine. Appendix A presents a discussion of the 
sampling routine which is being utilized and the probability distributions which are available 
for sampling purposes. 

In Figure 3.1 we start by setting the K counter equal to 1. This counter represents the 
number of frequency bands which are stepped through during the receive signal acquisition 
process. Details of this process are discussed later in this report. Next the counter i is set 
to zero. This counter represents the total number of samples (N) to be taken over the specified 
time period for the simulation run. In order that no information be lost due to the sampling 
interval the number of samples must satisfy 

N > 2T b W (8) 

.where T^ Signalling Interval (Sec.) 

W = Receiver Bandwidth (Hz) 
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Figure 3.1 INPUT SIGNAL GENERATOR 
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SIGNAL ACQUISITION MODEL 


4.1 


GENERAL 


For the MFSK system under consideration we must synchronize the signal in both time 
and frequency, where time synthroriization is defined as the correct timing between the 
transmitter and receiver at the beginning of each transmitted word and frequency 
synchronization is defined as the correct location at the receiver of the position of the 
carrier frequency in the spectrum. Figure 4.1 illustrates a basic MFSK receiver. The . 
receiver bases its decision upon an estimated Fourier power spectrum of the received 
input. It is assumed that an acquisition sequence of two alternating frequencies is 
transmitted prior to initiation of data transmission. 


4.2 


3 / 


FREQUENCY ACQUISITION AND TRACKING-' 


Equation (4) previously described the form of the received signal after carrier demodulation 
and will be repeated here 


X(t) = A cos (2^ f t + 0 } + n (t) 


( 4 ) 


4 / - 

. Using a discrete Fourier transform— to obtain the power spectrum of the received signal 
we calculate the Fourier coefficients 

2 2 1//2 


( 9 ) 


N-l 

= £ 

i=0 


X. Cos 

i 


(2 7T — ) 

V N ' 


k =0, 1, 2..., N/2 


00 ) 


3/ This discussion is drawn from JPL Space Programs Summary 37-52 and 37-57, Vol 

4/ Appendix B presents additional detail concerning the discrete Fourier Transform 
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Power Spectrum 



Figure 4.1 MSFK RECEIVER BLOCK DIAGRAM 












( 11 ) 


b k -g X,Sln(2-jt) l 

IT 

where X. = X(-jq-) are the sampled values of X(t). 


k =0, 1, 2. N/2 


The values of T & N are chosen with reference to the sampling theorem. If at least N 
samples are taken during the interval T (N = 2WT) where the process is low pass filtered 
to a bandwidth W and sampled at a rate 2W; then no information has been lost about the 
sampled waveform. 

Due to instabilities in the transmitter and receiver osci llators and to uncompensated doppler 
shifts, it is not possible to knew in advance the carrier frequency which will be received at 
the receiver. The actual received carrier frequency (f ) may be represented by f £ + X 
where X is a factor representing the difference between the actual and the nominal carrier 
frequencies. 

This initial frequency error is usually much larger than the bandwidth W of the Fourier 

spectrum process. However, an upper bound can be put on this initial error by a knowledge 

of the maximum anticipated doppler frequency shifts and by a knowledge of the long and 

short term oscillator drifts. This upper bound is say K W such that 

rr max 


-K W < f 

max ' c 


4 K W 
max 


If the local oscillator is stepped in frequency from f - ^ max ^ ^ c + ^ ma x^' ' n s * e P s ^ W 
and the spectrum, is computed in each band over a 2? second interval we can look for 
peaks in the band which exceed a threshold and use these peaks as a decision criteria for 
whether or not a signal is present. This is the procedure which is used in frequency 
acquisition. 

In order to track the frequency, we require an estimate of the frequency error. The 
estimator must be a function of the Fourier spectral components r, . One approximation 
to a maximum likelihood estimator of the true frequency used in this initial simulation 
is as follows. Other estimates are possible and should be tested . J. 

_ f +l - r "l 

- r ' TT 

+1 + - 1 

where ^ = maximum likelihood estimator at true Fourier frequency component 

r + .| = magnitude of the Fourier frequency component higher than the present estimate 

r .j = magnitude of the Fourier frequency component lower than the present estimate 

5/ "Frequency Tracking in An MFSK Receiver", by H. D. Chadwick, JPL Report 
#900-268, April 30, 1969. ; 


f 


This estimator can be used to design the frequency tracking loop. A general form of 
this loop is shown by Figure 4.2a. The estimator block output is defined as follows: 

t = E(^) + Nr = g (j) + N ? 

where: g (s) = E ( y ) is the nonlinear operation on the time error 

f produced by the estimator. 

= is the additional noise term. 

Figure 4.2a can then be rewritten as shown by Figure 4.2b. Substituting F(S ) = SH(s) 

Figure 4.2b becomes Figure 4.2c. It is seen that this is equivalent to the standard block 
diagram of a phase-lock loop with nonlinearity g( s ) replacing that of the phase detector. 

A second-order loop can be synthesized and because of the highly nonlinear nature of the 
function g( s ) the theoretical performance of the loop can be calculated by 

© The use of nonlinear approximation to the nonlinearity g( s ) with 

the standard linear phase-lock loop approximation. 

© non-linear phase-lock loop analysis 

The values of the variance of the frequency error, for the linear approximation and non-linear 
methods respectively are desired for comparative purposes. 

The tracking loop must convert the estimator into an estimate of frequency error A f. 

An actual MFSK receiver updates its estimate of the frequency error once every T seconds, 
requiring that the error remain small between these discrete sampling times. Three different 
discrete equivalent tracking loop filters for simulation purposes are shown by Figure 4.3. 

A loop simulation program can be implemented on a digital computer using Monte Carlo 
techniques to test each of the equivalent tracking loop filters. The objective is to obtain 
an interpolated error function (A f ) at a rate ten times faster than the calculation of the 
estimate. From the results, that version offering the best solution, will be implemented 
in the overall MFSK simulation program. 

4.3 TIME SYNCHRONIZATION-^ 

Again assuming a receiver bandwidth of W and sampling at or above the Nyquist rate as in 
equation (8) then if the receiver is not time synchronized with the transmitter, then the N 
samples will overlap two frequencies. Figure 4.4 illustrates the case. Note that as previously 
we are assuming an acquisition sequence of two alternating frequencies (f^ & is transmitted 
prior to data transmission. As seen in Figure 4.4, we have n samples from f^ and (N-n) 
samples from 


6/ Drawn from: "Time Synchronization In an MFSK Receiver", by H.D. Chadwick, 

J PL Space Programs Summary 37-48, Vol Ml. 
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Figure 4.4 TIME SYNCHRONIZATION ILLUSTRATION 


It has been shown— that a maximum likelihood estimator of the percentage time shift 
(E) is given by 

2-( r r-+ r 2) ; 

where r, & r_ are the Fourier coefficients from the F FT output. 

12 • 

Figure 4.5 gives a block diagram of a closed loop time synchronization technique which 
can be simulated by digital techniques. 


PREVIOUSLY DEVELOPED SIMULATION MODEL 


JPL has programmed the equations presented above. Figure 4.6 presents an overall flow- 
chart of the program and Figures 4.7 to 4.12 present one more level of breakdown for Figure 4.6. 
Work is underway to convert and modify this program for use on the Interdata machine. 

Appendix C presents detailed flowcharts of the pertinent sub-programs included in 
Figures 4.6 to 4. 12. 






































Figure 4.8 .. TIME ACQUISITION (See Figure 4.5) 
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Figure 4.9 DETECT 
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Figure 4.10 



FREQUENCY TRACKING AND INTEGRATION 
(See Figure 4.2C, Equivalent Phase Lock Loop) 






Compute 
S/N Ratio 


. I 



r 
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Figure 4.11 S/N ESTIMATE 
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Figure 4. 12 



TIME SYNCHRONIZATION TRACKING 
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APPENDIX A 
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DISCRETE FOURIER TRANSFORM CONSIDERATIONS 



A-l 


A 1/ 

DISCRETE FOURIER TRANSFORM CONSIDERATIONS - 


The Fourier transform has long been used for characterizing linear systems and for identifying 
the frequency components making up a continuous waveform. However, when the waveform 
is sampled, or the system is to be analyzed on a digital computer, it is the finite, discrete 
version of the fourier transform (DFT) that must be understood and used. The fast Fourier 
transform (FFT) is simply an efficient method for computing the DFT. The FFT can be used 
in place of the continuous Fourier transform only to the extent that the DFT could before, 
but with a substantial reduction in computer time. 


s(t) 
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Rgure A.l Discrete Venus Continuous Fourier Transform 


A.l/ Drawn from: "A Guided Tour of the Fast Fourier Transform" by G. Bergland, IEEE 
Spectrum, July 1969. 

"The Measurement of Power Spectra", by Blackman & Tukey, 

Dover Publications, 1958 

A-2 


Consulting Rgure A. 1 we illustrate the differences between the discrete and continuous 
transforms. Line (a) represents the continuous Fourier transform of S(t), a cosine wave. 

W(t) in line (b) represents a finite time window and its continuous transform through which 
we obtain samples of the signal. The portion of S(t) to be analyzed from the sampled data 
is shown in line (c) as the product of S(t) and W(t) in the time domain and as their 
corresponding convolution in the frequency domain. Lines (d) and (e) show the digital 
sampling operations and the resulting transforms. The continuous frequency domain 
function shown in line (e) can be made discrete if the time function is treated as one 
period of a periodic function (line (f)). Thus the DFT is simply a reversible mapping of 
N terms of a (k) into N terms of S\j). 

The three problems most often encountered in using the discrete Fourier transform appear 
to be aliasing, leakage, and the picket-fence effect. The term "aliasing" refers to the 
fact that high-frequency components of a time function can be impersonated by low frequencies 
if the sampling rate is too low. This uncertainty can be removed by demanding that the 
sampling rate be high enough for the highest frequency present to be sampled at least 
twice during each cycle. 

The "leakage" problem is associated with the effects of the rectangular data window shown 

in Rgure A.l line (b) resulting in a frequency domain function with an amplitude of the 
• 

— — form. The sidelobes of this function represent leakage into frequencies other than 
the true frequency. 

The "picket fence" effect is caused by the cumulative effect of the • Sm X spreading of each 

X 

Fourier coefficient. Thus instead of a line spectrum for each coefficient term we observe 
the phenomena shown in Figure A. 2 



1.0 

0.4 

0 


Power response 

rnmvm 



Rgure A. 2 The Response of the Discrete Fourier arid Transform Fourier Co- 
efficients Viewed as a Set of Bandpass Filters 
(Picket-Fence Effect). 
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APPENDIX B 


GENERALIZED RANDOM VARIABLE GENERATOR 



B-l 


PURPOSE 


The CURL subroutine may be used to generate random variables from several different 
standard distributions in addition to arbitrary, by inverse cumulative distribution curves 
read by the RDC UR subroutine. 


CURL 

The CURL routine is called by using '...CURL (ICV, Pi, @2...' in an arithmetic 
expression where ICV, Pi and P2 are described in detail in Table B.l 

The RNG routine is used by CURL to generate a uniform random number, N, between 0. 
and 1 . This random number is then used as the argument for an inverse cumulative 
probability distribution which may be formally evaluated or given in the linear curve 
format specified for the YLIN subroutine. The YLIN routine is used by CURL to evaluate 
the arbitrary curves. 

The distributions currently allowed in CURL are as follows: 


ICV 

>0 

0 

-1 

-2 

-3 

•4 

-5 

-6 

-7 

-8 


Distribution 

Arbitrary inverse cumulative probability linear curve [ICV] 

Constant 

Exponential 

Bernoulli 

Weibull 

Extreme val ue 

Normal 

Log normal 

Uniform 

Beta 
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CURL PROBABILITY LAWS AND DISTRIBUTIONS 
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CURL PROBABILITY LAWS AND DISTRIBUTIONS (Cont) 



Beta ' a > 1 b>l .. _ J^a + b}. ®, r n * , See page 945-26 . 5 . 22 of 

T(a)Y{b) * u < t<: 1 NBS Handbook of Math 

= 0 elsewhere Functions 
















RDCUR 


The RDCUR routine may be used to read in any inverse cumulative probability linear 
curves which are to be used by CURL in the run. When RDCUR is called, the following 
data are read from unit 5. 


Card Col umns 


1 1-10 


2 1-10 

3 6-10 

col. fields 


Repeat card 3 as needed. 

Repeat cards 2 and 3 as needed. 


Description 

NCUR, number of arbitrary inverse cumulative 
density curves (^0). If NCUR-^0, the follow 
ing cards are read for each curve: 

NPTS, the number of points in this linear 
curve 

The coordinates of the curve points consisting 
of the probability, followed by the value of 
5, in non-decreasing order of probability. 

(3 points per card). 


The above input data if NCUR is greater than zero must be stored in the following 
COMMON variable locations which are accessible to the CURL routine. 

NCUR. ...... holds the number of arbitrary curves 

CURL ( )... dimensioned sufficiently to hold the tables 

LOCCUR ( ). dimensioned at least NCUR. 

Subroutine used by CURL 
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RNG 

YLIN 

XNORD 


APPENDIX C 


SIGNAL EXTRACTION FLOW CHARTS 
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APPENDIX G- SIGNAL EXTRACTION FLOW CHARTS 


The flow charts contained in this appendix were developed by analysis of the program 
listings for the JPL simulation of the MFSK low bit rate receiver. Neither program 
documentation or program flow charts were available other than the listings. 

The contract requires that the JPL simulation be utilized on the Interdata machine 
as a portion of the baseline simulation being developed. Thus, analysis of the 
program listings and development of detailed flow charts is a necessary first step in 
understanding what was simulated and in converting the JPL programs written in both 
FORTRAN and machine language to the language used on the Interdata machine. 

The mathematical description of the time and frequency synchronization operations 
(described in the main body of the report) were also utilized in developing an under- 
standing of JPL program listings. 
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CALL FREQ (I FREQ) 
CALL NIXIE (IFREQ) 



CALL TIME AC (TCI) 


ERASE 

MICKEY(MOUSE) 

PAUSE 

ERASE 

ALPHA 350, 225, 100,0 
TIMEACQU. MODE XX 


EOMODE 030010 
1777 0 (ISTART+ 2400) 

' ^ I START 

17770 (ISTART+ 1001 

IFIN 

(IF1N+ 768) O 1777 
— (HALF 


IFLAG 

JFLAG 

NCOUNT 

NERROR 

ISIG1 

I SI G 2 

FRINT 

FERR 

SUMS 

SUMN 


© 


Figure C-l (Conflnued) 
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’Figure C.l (ConHnued) 
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(I START + 3000) O 1777 
ISTART 

(ISTART+ 1001) O 1777— ^ 1FIN 
(IFIN + 768) n 1777 — ► I HALF 






Figure C.l (Continued) 
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INITIAL FREQ ACQU. SYNCHED LOC. OSC. FREQ. IS 

RETURNED AS INFREQ. IF LOCK IS OBTAINED, LOCK = 1 

IF NOT, LOCK =0 



Figure C.2 FREQAC (ALPHA, IF RE o’; LOCK) 

























0-*-K 

0— MRAG 
EOM 020020 
POT = 02004101 



\ 

i—i 

) — ; 

, ! 



STOP 




NT^\ 



YES 


CALL VVFC 
IDATA (ISTART)L 
IQUT, ISIN, 9j 
ISCRATCH, 1) 



Figure C.3 TIMEAC (TCI) INITIAL TIME SYNCH AND TABLE SETUP 
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Figure C.3 (Continued) 
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DETERMINES LARGEST SPECTRAL COMPONENT IN THE RECEIVED SPECTRUM 

ISIG1 !S THE LATEST SIGNAL 
ISIG2 IS THE PREVIOUS SIGNAL 



Figure C. 4 DETECT (ISIGI, ISIG2) 







determines largest spectral component in the received spectrum 

ISIG1 IS THE LATEST SIGNAL 
ISIG2 IS THE PREVIOUS SIGNAL 



Figure C.4 (Continued) DETECT (ISIGI, ISIG2) 
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UPDATE FREQ. TRACK LOOP. INTEGRATION IS 
PERFORMED BY THE INTEGRATE SUBROUTINE 



Figure C.5 FREQTR (IS1GI, IFREQ, FRINT, FERR) 
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(SUMN * (NCOUNT-1) + AVNO)/ NCOUNT — - SUMN 



Figure C.6 SNEST (IWORD, SUMS, SUMN, NCOUNt, ISNR) 
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TIME SYNC TRACKING 



(IOUT (ISIG 1)+ IOUT (ISIG 1 + 256) 2 ) * IFILTER(ISIG 1) 

RSQ 1 

SAME F OR ISI G2 ► RSQ2 

VRSQI R1 

VRSQ2 * R2 

256 * (R2-R1)/(R1 + R2) . 

► EPS 

(START + TC2 * EPS START 

START (START 


l 



Figure C. 7 TIMETR (TC2, ISIG 1, ISIG2) 

• ' • * 
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INTEGRATION FOR THE FREQUENCY TRACKING LOOP 
RATE CONTROLLED BY INTERRUPT 215 (SET TO 1/SEC) 
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Figure C.9 BLIP 
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TOTH P E A ADC INra C RUPT ULTIPLEX CHANNEL0 ' IDA ™ IS CONNECTED 



IDATA,2 

— ^ IDATA+ 1024,2 

- 

i 

i 



EOM 030064 
PIN IDATA.2 



©— ' 



i 



IPOINT n 1717 
— IPOINT 
300 S — «- 20 IS 
SET SWITCH OFF 

,, YES^ 



^START^) 



EOM 030066 
EOM 030070 
POT = 0 

100S — c- 207 (SET INTERRUPT TO 200S) 
301S — ► 2015 (SET SWITCH ON) 

1023 — >- (POINT 

512 SHIFTED — •- IVERT 

EOM 030031 

POT IVERT 

EOM 030003 

EOM 030005 


^RETURN^ 



STORE A, B, X 
IPOINT+ 1 — *• 
SHIFT 12 LEFT 


IPOINT 


Figure C. 10 SETUP 
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Figure C.13 MICKEY (MM) 
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